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Abstract 

New lower bounds for the binding energy of a quantum-mechanical system of 
interacting particles are presented. The new bounds are expressed in terms of 
two-particle quantities and improve the conventional bounds of the Hall-Post 
type. They are constructed by considering not only the energy in the two- 
particle system, but also the structure of the pair wave function. We apply 
the formal results to various numerical examples, and show that in some cases 
dramatic improvement over the existing bounds is reached. 
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I. INTRODUCTION 



One of the central problems in quantum many-body physics is to find the energy of a 
system of A particles interacting with given two-body potentials. Variational techniques 
yield an upper bound to the exact energy. The determination of a strict lower bound can 
then provide a natural and useful complement. 

In recent years there has been renewed interest in deriving such lower bounds, mostly 
in connection with quark models in hadron spectroscopy or the limits for Borromean 
phenomena in loosely bound systems @-||]. Earlier uses of lower bounds were focussed on 
thermodynamical considerations ||], or the problem of stability in self-gravitating systems 

01- 

Up to now all lower bounds are based on the Hall-Post decomposition P,pH[] of the 
hamiltonian into two-body clusters, and subsequent application of the variational principle 
in two-body space. Several variants and extensions have been proposed, e.g. an optimal 
decomposition in case of three |T1] or four |12j unequal masses. The case of identical fermions 



was recently studied in [13 



In its most useful form, a lower bound of the Hall-Post type is expressed in terms of 
the ground-state energy in a two-body system. Finding this two-body energy usually goes 
together with determining the wave function of the ground-state pair. We show in this 
paper that by using the structure of the pair wave function one can always improve the 
Hall-Post bound, and that the improvement is sometimes spectacular. Our results apply to 
boson and fermion systems, both with and without the presence of an external potential, 
and irrespective of the local or nonlocal character of the two-body potentials. We do restrict 
ourselves to systems of identical particles, leaving a study of the unequal mass case for future 
work. 

Loosely speaking, the Hall-Post bound implies that the ground state for A identical 
particles is a superposition of \A(A — 1) pairs in the lowest-energy state of a modified 
hamiltonian. However, because of the correlated structure of such a pair, it is in general 
not possible to reach a pair occupation of |^4(^4 — 1), except in the case of noninteracting 
bosons. For any given pair, there exists a maximal value for the occupancy that can be 
reached in an A-body state. The remainder of the \A(A — 1) pairs must then necessarily 
be in a state with higher energy. This allows the construction of a new lower bound, whose 
value is increased with respect to the Hall-Post case. 

The present paper is organized as follows. In Section |J we first establish notations, 
introduce cluster decompositions of the many-body hamiltonian, and rederive the Hall-Post 
inequalities. Next we point out in Section [Hj how the Hall-Post bound can be improved 
in the case of fixed-center bosonic systems. The equivalent case for fermions is treated 
in Section |V[ The modifications that enter when considering self-bound (translationally 
invariant) systems are discussed in Section |V[ Finally, we apply the formal results to var- 
ious numerical examples, which are collected in Section [VI]. Section |VI]] contains a global 
summary and points out some remaining problems. 



II. GENERAL REMARKS ON CLUSTER DECOMPOSITIONS OF THE 
HAMILTONIAN FOR A SYSTEM OF IDENTICAL PARTICLES 
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A. Notational conventions 



We restrict ourselves to systems of identical particles and follow the notational conven- 
tions of ||14|| . In particular, single-particle coordinates are generically written as X\ 7 X2 • • ., 
and should be regarded as spatial coordinates in d dimensions. Spin or isospin degrees of 
freedom are not explicitly mentioned, but may be assumed to be included in the Xi depen- 
dence of the wave functions and operators. 

For a system of A particles we consider a hamiltonian //^[a] which is a sum of one-body 
and two-body operators with variable relative weight, 

H A [a)(x 1 ,...,x A ) = a ^ t(x h ) + v(x h ,x i2 ). (1) 

h=l,...,A ii<i2=l,..., J 4 

Using a combinatorial identity we can rewrite H A \a = 1], for any iV with 2 < N < A, as 

H A [l](x 1 ,...,x A ) = ( A \ 2 \ E H N [——^]{x h ,...,x iN ). (2) 

J h<...<i N =l,...,A 

The hamiltonian ifjy[^5f ] f° r clusters of N particles is related to the original hamiltonian 
H A [1] for A particles, but has different relative weights for its one-body and two-body 
components (or, equivalently: a different coupling strength of the two-body interaction). 
We can now make the weights equal again by absorbing the difference in an additional 
one-body term. As shown in Section [II B , this leads to sumrules for the energy which 
are generalizations of the well-known Koltun sumrule [IB]. Alternatively, we can keep the 
different coupling strength in the iV-particle hamiltonian. This will lead to another set of 
sumrules, derived in Section HIC , from which the Hall-Post inequalities immediately follow. 

One of the ingredients in these sumrules are the spectroscopic factors related to the 
removal of particles from the A-particle system in its ground state. These are defined as 
follows)]. 

The A-particle ground state \I/o(yi)[l] of ^Ta[1] can always be expanded in terms of the 
complete orthonormal set of iV-particle eigenstates ^(jy) [a] of H N [ot\, 

— 1 /2 

*o(^)[l](^i, • • • ,x A ) = (n) E ^iHf 1 !- • • • ,x N )ip u ( N )[a}(x N+lj . . . ,x A ). (3) 

u(N) 

The expansion coefficients ip u (N)[a] are the overlap functions between the ^(jvJai] and 

*0(X)[1], 

4> u (N)[(y}(x N+1 , . . . ,x A ) = (4) 

1/2 



j dx l ...dx N % {N) [a](xi,...,x N )V ( A) [l}(x 1 ,...,x A ) 
and their normalization yields the corresponding spectroscopic factor S u ^)[a], 



The definitions in this Section may look tedious or superfluous, but will allow in later Sections 
a unified treatment for fixed-center and self-bound systems. 
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S v (N)[a] = J dx N+ i . . .dx A \ipv(N)[a](x N+1 , . . . ,x A )\ 2 , (5) 

that is, the squared amplitude for removing A — N particles from \E f (A)[l] an d ending up in 
^u(N) [«]■ 

Obviously, the spectroscopic factors are positive, and the completeness of the set ^(at^o;] 
implies the sumrule 

S HN)[®} = (n)- (6) 

u(N) 

Note that the S^ao [«] can also be interpreted as the occupation number of the state 
^(jv) [at] in the ground state ^(A)^]) since Eq.(||) can be rewritten as 

S„(N)[oi] = J dxi- . .dx N dx[ . . .dx' N ^l {N) [a](xi } . . . , x N )p^ N \x u . . . } x N ]x' l} . . -,x' N ) 

^(aoHK,...,^), (7) 

in terms of the iV-body density matrix 

p ( - N \x A -N+i, ■ ■ ■ ,x A ;x' A _ N+1 , . . . ,x' A ) = J dxx...dxA-N (8) 
^(ajWC^Ij ■ ■ ■ ,x A - N ,x' A _ N+1 , . . .,x' A )^ ( A ) [l](a?i, • • .,x A -n,xa-n+i, ■ ■ -,x A ). 

The above equations hold for fixed-center systems, in which the one-body part of the 
hamiltonian contains an external potential. For self-bound systems the one-body part is 
purely kinetic, and the two-body part depends on relative coordinates only, i.e. v(xi,x 2 ) = 
v(xi — £2). Since the intrinsic eigenstates of the hamiltonian are translationally invariant, 
some modifications [14]] are needed. 

The intrinsic kinetic energy is obtained by subtracting the center-of-mass from the total 
kinetic energy, 

Am h=l,...,A A/im h=l,...,A A/im h<i 2 =l,...,A 

and is seen to behave as a two-body operator with an A-dependent coupling strength. 
Defining the intrinsic hamiltonian with variable coupling strength as 

H' A [a] = af' A + + v(x h -x i2 ), (10) 

i 1 <i 2 =l,...,A 

the decomposition (0) gets modified to 

1 - N 

H' a [1](xi,...,xa) = , A2 x E H' N [—](x hl ...,x iN ). (11) 

( JV-2 J ii<...<tjv=l,— i-A 

The A-particle ground-state wave function can be expanded in a complete orthonormal 
set of intrinsic iV-particle eigenstates |b| , the self-bound equivalent of Eq. (^) being 
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— 1 /2 

u(N) 



xip„( N )[a](x N+ i - R N , ...,x A - R N ), (12) 



where Rn = j? Z)£Li x u the overlap functions are defined as 



i>u(N)[a]{ x N+i, ■ ■ ■ ,x A ) = ( n) J dx x . . . dx N 5(R N )%( N )[a\(xi, . . .,x N ) 

x^f 0{A) [l}(x 1 ,...,x A ), (13) 

whereas Eqs.(|]||) remain unchanged. 

The spectroscopic factor Sum, as defined by Eq.fl^), can again be viewed as the occu- 
pation number of the intrinsic iV-particle state ^ u (n) in the A-particle ground state, but 
this occupation number now also includes a summation over all possible states of iV-particle 
center-of-mass motion, i.e. using 

we can express the spectroscopic factor as 

Sv(N)[a] = J dx 1 . . .dx N dx[ . . . dx' N 5(R N - R' N )^l^[a](xi, . . .,x N ) (15) 
p (Ar) (xi, . . .,x N ;x[, . . .,x' N )^„ {N) [a](x[, . . .,x' N ). 
The iV-body density matrix for a self-bound system, defined as 

p ( - N \x A ^ N+1 ,...,x A ;x' A _ N+1 ,...,x' A ) = J dx x ...dx A ^ N 5{R A - N ) (16) 

^(^[^(^l) • • • > X A-N: X A _ N+1 , . . . ,X^)\l/o(A) [l](a?l, • • • , ^A-AT, X A -N+U ■ ■ ■ ? x a), 

is the proper extension of the one-body density matrix for translationally invariant systems 



m 



14,15 



B. Generalized Koltun sumrules 

For fixed-center systems the decomposition (0) can be reshuffled so as to have the same 
hamiltonian in A-particle and iV-particle space, at the expense of an additional one-body 
term, 

H A [l](x l ,...,x A ) = , 1 v H N [l](x il ,...,x iN ) - — N t(x h ). (17) 



<...<i N =l,...,A ii=l,...,A 



The expectation value of ([17]) in the A-particle ground state, expanded according to 
Eq.(H), can be worked out using Eqs.(f|,[D and reads 

£o ( A)[l] = aTW^T ' To< ' 4) " T1~T\ E (S-otM - J3o(.4)[l])S„ (K) [l] \ . (18) 



(j--0 
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Here the ground-state energy has been expressed in terms of the expectation value T ^ of 
the one-body field and the first energy-weighted moment of the distribution of (A — N)- 
particle removal strength (or mean removal energy). This is a generalization of the familiar 
Koltun sumrule |Tj| for N = A — 1, 



Eo{A)[l] 



1 



To(A)- E {EuiA-Dm-Eo^S^A^ll}}, 

u(A-l) 



(19) 



where the distribution of single-particle removal strength is experimentally accessible 
through single-particle (SP) knock-out reactions ||T7| , [l8| . 

For self-bound systems we can [in analogy to Eq.([H])] restore in Eq.flTTD the original 
hamiltonian at the expense of an additional (intrinsic) kinetic term, 



H' A [l](x u ...,x A ) 



y n - 2 j i\ 



A — N - 

E H' N [l](x h ,...,x iN ) t^T' a . 



<...<i N =l,...,A 



N 



(20) 



The ground-state expectation value of Eq. (|20"|) , evaluated by means of Eqs.(fi|,|l3D, then 
reads 



E o(A) [1] 



N 



1 



52(E v[N) [l}-E 0{A) [l})S uiN) [l} . (21) 



u(N) 



A + N-l 
The result for the case of N = A — 1 

^0(A)[1] = - | A _ 1 T Q ( A) - Yj (^(A-l)[l] - E 0(A) [1])*S'^(A-1) [1] 



(22) 



is the Koltun sumrule for self-bound systems (including the correct recoil factor for the 
kinetic energy [14, HI). 



C. Hall-Post inequalities 

The expectation value of Eq.(Q), with the modified coupling strength in the iV-particle 
hamiltonian, can be evaluated in the same way as Eq . fll8|) , 

1 ^ N — 1 N — 1 

[n-2 ) u{N) /ll /ll 

This expression forms the basis for deriving the lower bounds considered in the next Sec- 
tion. At the simplest level, combining Eq.(^3|) with the trivial inequality resulting from the 
spectroscopic sumrule @, 

0<^ w [^j]<(£), (24) 
immediately leads to inequalities of the Hall-Post type 
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m^rf^ 1 ^- (25) 



For self-bound systems the A-particle ground-state energy can be likewise evaluated, from 
the expectation value of Eq. (|TTD , in terms of intrinsic iV-particle energies and occupation 
numbers, and reads 

Eo(A)[l] = r^jY E E V(N) [^}S V(N) £}. (26) 

u(N) A A 



(n-1) 

The Hall-post inequalities 



A(A - 1) N 
E o(A)[l\ > jv(iV- 1) 0(JV) ^' i2l) 



then follow immediately from Eq.( 



III. IMPROVED LOWER BOUNDS FOR FIXED-CENTER BOSONIC SYSTEMS 

A. Derivation of new lower bounds 

From now on we concentrate on the most relevant case N = 2, and try to derive lower 
bounds for the A-particle ground-state energy in terms of two-body quantities. At present, 
such a lower bound is given by the Hall- Post inequality (|25|) , 

E oi A)[l]> (^Eowij^}, (28) 

which is usually derived by applying a variational principle in two-body space ||. As was 
shown in Section [11 Q the Hall-Post inequality is actually an equality when expressed in 



terms of two-body energies and occupations, 

E o(A)[l] = E^Ij^lWj^y]- ( 29 ) 

In order to simplify notations we drop the coupling strength dependence in this Section, 
as it will be fixed to a = 1 for A-particle and a = -j—r for two-particle quantities. It is 
clear from Eqs.( |28"| , |2"9| ) that in the traditional lower bound the distribution of pair strength 
S v {2) is approximated by concentrating all strength in the two-body ground state, i.e. by 
the distribution 

S v (2) ~ <5o^( 2 ), (30) 

which exhausts the sumrule fl|). 

This assumption can only be realistic for weakly correlated systems. For a non-interacting 
Bose system [v = in Eq.(|l])] it holds exactly. The uncorrelated eigenstates are product-type 
wave functions, 
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^ ) {x 1 ,...,x A ) = II Xfa), %$(xi,x 2 ) = xMx{x 2 ), (31) 

i=l,...,A 

where x( x ) 1S the SP eigenstate of the one-body hamiltonian t(x) corresponding to the lowest 
energy 6q. As a consequence, SqSK = A(A — l)/2, and the Hall-Post lower bound coincides 
with the exact result, 

E o(A) = Ae o = ( a ) 7 _ x 2 eo = ( 2 )^o(2)- ( 32 ) 

For strongly correlated systems we expect that the occupation Sqm) can deviate substantially 
from A(A — l)/2, and the Hall-Post lower bound will be far from the exact result. 

In order to improve on this situation we need to take into account correlations in the 
structure of the many-boson eigenstates. As it is more convenient to work in second quan- 
tization, we define </?L 2 ) as the creation operator for the two-body eigenstate ^f u (2), i.e. 

I *,(2)> = V«/(2) I)- (33) 
The occupation if>„(2), as defined by Eq.(|j) or Eq.(0), is then simply written as 

Sv{2) = (*0(A) I rf(2)W(2) I *0(A))- (34) 

We now consider the upper bound S^ 1 ^ for S„ip) , 

Sffi = max { (9 A I rf (2) ^ (2) I * A ) } , (35) 

where the maximum is taken with respect to all normalized A-boson eigenstates. This 
upper bound only requires knowledge of the two-body state ^(2); its explicit construction 
is pointed out in Section [111 B . 

Since S^f? is better than the trivial upper bound fl2"4|) , the following inequality holds, 

Sv(2)<S^<( A 2 ). (36) 

In combination with Eq. (|29|) this results in a sequence of new lower bounds £ M , with \i = 
0, 1, . . . , for the ground-state energy, 

f- 1 max f ^-1 1 

Eq( A ) >£v = Y1 S ^(2) E u(2) + < ( a J - S ™2) ( ( 37 ) 

Here the two-body states are assumed to be ordered according to increasing energy, E / 2 ) < 
-Ei (2) < — 

The optimal lower bound in this sequence is given by the largest value Cp, = max /t (£ /i ), 
where Jx is determined by 

u=0 u=0 
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The inequalities in Eq.fl3~7|) constitute the principal result in this paper, and several 
remarks are in order: 

(a) The conventional Hall- Post bound of Eq.(|28|) coincides with £ . 

(b) Ci is always a better bound than £ , since 

A - Co = {( i ) - S^} [E m - E m ] > 0. (39) 
In most of the cases we studied £1 is the optimal bound C^. 

(c) If only a finite number rib of discrete levels is present in the two-body spectrum, then 
Eq.(|3~7|) holds for < /i < rib, where E nb can be taken equal to zero. 

(d) Eq.(|3"TD holds without any symmetries of the underlying hamiltonian. If these are present, 
they can of course be used to refine the lower bound, e.g. an energy level E v m) appearing 
in Eq.(^) can correspond to a d^-fold degenerate multiplet with eigenstates <£>„(2) jA t> A* = 
1, . . . , d u . If the quantum numbers of the ground state ^o(A) are known, it will usually be 
possible to determine the maximal joint occupation number of the multiplet, 

= max | (* A | £ iplw^w* I *a> } , (40) 

where the variation is made over all A-particle states with the same symmetry properties 
as ^o(A)- Using P^f m the evaluation of Eq. (|37|) will result in a better bound, because 
Pu{2) — S™(tyn- ^n example of this will be given in Section [VI]. 

(e) Finally, doing better than would require e.g. an optimalization of the simultaneous 
occupation for the two-body ground and first excited state, involving a determination of 

max(^ A | #0(2)^1,(2)^0(2) + #1(2)^1(2)^1(2) I *a)- (41) 

This appears to be a far more complicated problem than the determination of S^y 



B. Construction of maximal pair occupancies 



Introducing the set of its natural orbitals Xa{x), a general two-boson state \I/(2) can be 
written as (see Appendix |A 1| ) , 

V(2){X 1 ,X 2 ) =^2x a Xa(x 1 )Xa(x2), (42) 
a 

where the x a are real and positive, and J2a x a = 1- 111 secon d quantized form this reads 

\^(2)} = -^T,^(ci) 2 \) = V \ 2) \), (43) 

where is the creation operator for the one-boson state Xa(x). 

The maximal pair occupation S^ x of ^(2), defined according to Eq . (|3~5D , is equal to the 
largest eigenvalue of the following hermitian eigenvalue problem in A-boson space, 

¥>(2)¥>(2) I = A I (44) 
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The corresponding eigenvalue problem for fermions was recently solved by Pan et al. 
2"0fl , in the context of a generalized pairing problem. The method in [20 ] , which involves an 



infinite-dimensional algebra, can be easily adapted to bosonic systems. Here we only state 
the final result; for completeness the derivation is given in Appendix 0. 

Let s = for A even, s — 1 for A odd, and a = (A — s)/2. The largest eigenvalue of 
</>( 2 )¥>(2) 1S among the eigenvalues A in the maximally paired subspace (see Appendix ||). 
These are given by 

A = 1 + 2*4 + 4 J2 — > ( 45 ) 

fc=i,...,o-i y k 

in terms of the solutions of a set of a — 1 nonlinear equations in a — 1 variables (yi, . . . , y a -i), 
1 ^ x 2 a {l + 2s5 a p) 1 ^ 1 

tLh 2 — = — + E > fc = l,...,a-l. 46 

4^ l-3/fcx2 y fc i=i,..t^i-MkVk-yi 



The index /3 which appears for odd A is arbitrary (see Appendix |B|) . 

The system of equations ( f46|) allows to determine the maximal occupancy for general 
A. For small values of A it can be rewritten in a much simpler form: 



(A = 3) — > A = 1 + 2xJ, 

V A - 1 - 4^2 

The maximal occupation of a pair state depends on the structure of ^(2), i-e. on the 
distribution of the x a . The two extreme cases are: (a) uncorrelated, and (b) "maximally 
correlated". In the uncorrelated limit (a) only one of the x a is non-zero, resulting in S^ x = 
A(A — l)/2. The latter limit (b) has a flat distribution, i.e. assuming that there are Q SP 
states, then x a = This case, corresponding to a schematic boson pairing force in a 
single degenerate shell, can be treated analytically since the algebra reduces to SU(2). The 
resulting maximal eigenvalue is S"^ ax = a(l + ) • 

In summary, the maximal eigenvalue in A-boson space obeys 

a<S^<( A 2 ), (47) 

where the upper limit corresponds to an uncorrelated wave function, and the lower limit is 
reached as the large-f2 limit of a maximally correlated wave function. 



IV. MODIFICATIONS FOR FIXED-CENTER FERMIONIC SYSTEMS 



The basic inequalities (|37|) derived in the previous Section still hold for fermions. How- 
ever, the boundaries for the allowed pair occupation numbers are completely different. 

The natural orbitals for a two-fermion state come in associated pairs (Xa(%), Xa( x )), an d 
a general two-fermion state \&(2) can be written as (see Appendix |A2|) , 
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^(2)Ol,X 2 ) = X a \Xa{xl)X&(x 2 ) - Xa(xi)Xa(x 2 )}, (48) 

where the sum runs over distinct pairs. The x a are real, positive for a > 0, and with 
Xa = —x a , normalized as J2a x a = 1- m second quantized form this reads 

|* (2) ) = V2E^4l) = ^ 2) l)- (49) 

The maximal pair occupation 5*^ x of \l/(2) is again equal to the largest eigenvalue A of 
an eigenvalue problem in A-fermion space, 

¥>(2)¥>(2) I *a) = A | (50) 

For fermion systems we can use directly the results in [2C|. Let s = for A even and 
s = 1 for A odd; and a = (A — s)/2. We have to solve the set of a — 1 non-linear equations 
in a — 1 variables (yi, . . . , y a -i), 

_i £ ^a-,M = 1+ E _^ ii = li ... !0 _ L (51) 

^ a>0 1 ^ 2=l,...,a-l;Z^fc ^ 

Then the relevant eigenvalues of (|50| ) corresponding to the subspace of maximally paired 
states are given by 

A = 1 - 2sx 2 p - 4 J2 — ( 52 ) 

fc=i,..., -i y k 

The simplified secular equations for small A read, 

1 - 2x% 





3) - 


-> A = : 


(A = 


4) - 




(A = 


5) - 


«>o 



A - 1 + 4x2 

a£(i-M 

A - 1 + 4x2 + 2x| 



1. 



It is again instructive to consider the two limiting cases in the structure of ^( 2 )- In the 
uncorrelated case only one of the coefficients x a is non-zero. This corresponds to a two-body 
Slater determinant, and 5"^ ax = 1. In the maximally correlated case, where all coefficients 
are equal, the coefficients become | x a \— -7=, if there are Q SP states. This is equivalent 
to the well-known problem of a schematic fermion pairing force in a single degenerate shell, 
with S$ x = a(l- 2(Q ~ n 1+s) ). 

In the large-f2 limit we then find that the maximal pair occupation in A-fermion space 
obeys 

1 < S^ x < a. (53) 

Note the different role of correlations for boson and fermion systems: for boson systems 
correlations decrease the maximal pair occupation compared to the uncorrelated case, for 
fermion systems they increase it. 
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From these bounds on S^f* it follows that the Hall-Post lower bound ( p8|) will never 
be satisfactory for fermion systems. Even without knowing the structure of the two-body 
eigenstates, it can be replaced by the better bound 

E 0{A) >aJ2 E v(2), (54) 

v{2) 

where v = 0, . . . , 2 (a — 1 + s). Unfortunately, this bound does not in general become the 
exact result in the limit of noninteracting fermions. For a noninteracting fermionic system 
the new bound (|3"T| ) would yield = J2 V (2) Eu(2), where v — 0, . . . , \A(A — 1) — 1; hence the 
new bound, though better, would still not be exact, the reason being that the \A(A — 1) 
two-particle energies Eij = e« + £j made with the A lowest SP energies ej, are not necessarily 
the lowest two-particle energies. 



V. MODIFICATIONS FOR SELF-BOUND SYSTEMS 

In analogy to the treatment in the previous Sections we try to derive lower bounds for the 
intrinsic A-body ground-state energy in terms of the two-body wave functions and energies 
of relative motion, by considering Eq.fl2"E|) with N = 2, 

E (A)[1] = E E ^)& S »(A ( 55 ) 

u(2) A A 

We again drop the dependence on the coupling strength, since it will remain fixed at a = 1 
for the A-body, and at a = -f for the two-body quantities. 

Apart from the different couplings there are no differences with the fixed-center case, and 
the basic set of inequalities Eqs.flSTD is still valid. The novel complication lies in deriving an 
upper bound Sffi* f° r the P & i r occupation of a relative pair wave function \l/(2) m an intrinsic 
A-particle wave function ^(a)- Mathematically this boils down to finding the absolute 
maximum of the pair occupation in Eq.(|5|) or Eq. (|15D , 

S(2) = ( n ) J dxi dx 2 dx[ dx 2 5(R2)5(R' 2 )^ 2 - ) (xi, x 2 )^ (2)(x[, x' 2 ) 

fd*,... X 2 , l3 , . . . , 4, -3, ■ ■ ■ , X. 4 ) (56) 

by varying ^ (a) in the space of translationally invariant wave functions of the correct 
(anti) symmetry. 

The problem of spectroscopic factors and occupation numbers in self-bound systems is a 
difficult one (see and we did not succeed in finding a general solution to this problem. 
The case A = 3 is tractable, however, since it can be transformed into an eigenvalue equation 
in SP coordinate space. We neglect (iso)spin degrees of freedom and only consider here cases 
where the spatial part of the wave function is totally symmetric (rj = 1) or antisymmetric 
(77 = — 1) , though the results can probably be extended to cases of mixed spatial symmetry 
0. 

Introducing Jacobi coordinates a = x\ — 22 and b = x% — R 2 , a general three-body wave 
function ^/^) can be written as 
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a 3a 6, 



a 3a 6, 



*(3)(zi, x 2 , x 3 ) = * (3) (a, 6) = /(a, 6) + 77/(6 + - , — - -) + f(b - -, -— - -), (57) 



4 



in terms of a function f(a, b) = r]f(—a, b). The pair wave function is simply 

*( 2 )(ari, ar 2 ) = 0(a) = r)g{-a). 
In terms of these quantities, the pair occupation fl56|) for A = 3 is rewritten as 



(58) 



S(2) =3 / dbdada' f*(a! ,b)g*(a) 



<7(a')%) (a, 6) + W (| + 6)%) (a, ^ - ^) 



, a' . N ~ , 3a' 6 . 

+^(__ +&) Vl/ ( 3 ) ( G) -_--) 



(59) 



and the maximum must be taken with respect to all f(a, b) having a fixed normalization 

dxidx2dx3d(Rs) | \E f ( 3 )(xi, X2, £3) 1^=3 y dad6 /* (a, 6)^(3) (a, 6). (60) 
Performing the variation leads to the secular equation 

A§ ( 3)(a\6) = J dag* (a) g(a')V i3) (a, b) + W (| + 6)* (3) (a, ^ - ^) 

(61) 



.a' ~ 3a' 6 , 

+ ^(__ + & )^ (3) ( a) ____ y 



Introducing the overlap function G{b) = J dag*(a)^(3)(a,b) [see Eq. fljjD ], the secular equa- 
tion is transformed to a SP eigenvalue equation of a hermitian non-local operator, 



(A - l)G{x) = 2 V J dx' 



9*(^ + x')g(j + x) 



G(x'), 



(62) 



which can easily be solved numerically. 

Although we cannot yet determine the maximal pair occupation for general A, we can still 
find a bound for -Eo(A) i n terms of two-particle quantities that is better than the traditional 
Hall- Post bound Eq . $27\) for N = 2. This is done simply by replacing in the Hall-Post lower 
bound (0) for iV = 3, 



E 0(A) [a = 1] > 



MA -i) 

6 



E, 



0(3)1 



A 1 



(63) 



the three-body ground-state energy £"0(3) [4] by an improved lower bound obtained by solving 
Eq.©. 



VI. NUMERICAL EXAMPLES 
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A. Trapped boson system with pairing forces 



As a first example we consider a system of spinless bosons trapped in a (three- 
dimensional) harmonic-oscillator well, and interacting with a general monopole pairing force 



see, e.g., Dukelsky et al. [21]). The hamiltonian reads 



H = J2e N J2J mm c Nlm + gPlP, (64) 

N Im 

where N = 0, 1, . . . is the harmonic oscillator (HO) quantum number and I the orbital 
angular momentum. For convenience we remove the zero-point energy from the single- 
particle spectrum and put too — 1, i.e. we take = N . 
The pair operator in Eq.flBjl) is 



p ] = E wn E (4 ■ 4) = E wn E(-i)' V2IT1 [4 ® 4] ° • (65) 



N I N 



For the purely schematic pairing force, with constant wn = w, the system is exactly solvable 
for any finite number of HO levels, as was demonstrated by Richardson p2| . However, 
the schematic force has some unrealistic features due to its implicit dependence on the 



degeneracy = (N + 1)(N + 2)/2 of the HO shells [21]. The interaction between boson 
pairs in N and N' levels is proportional to y/D^D^i. For attractive pairing, e.g., this leads 
to occupations of the higher levels far exceeding those of the lower ones. This can be cured by 
taking wn = which is the pairing force we will consider in our numerical examples. 

In order to have the same notation as in Section TlfB, we can go over to the natural 
basis for P* by defining new SP states, 

4m = -^i"^ 1 -^)) (4 m + sgn(m)4_ m ) , if m + 0, 

4o = 4o; (66) 

in terms of which P^ = Y^Nim w iv(4m) 2 - 

The construction in Eq.([Tj]) of a lower bound for the A-boson system requires first 
to solve the two-boson problem with the same hamiltonian (|4]) but modified SP energies 
e N = ^Ej e N, and we briefly discuss its solution. The two-body eigenstates of the collective 
pairing type (involving all harmonic oscillator shells) can be written as 

I^(2)) = 4kE^™(4J 2 |), (67) 

V ^ Nlm 

and have eigenenergies E( 2 ) which are solutions of 

-1 = E WlDN , • (68) 

The corresponding wave function is x^im ~ wn/(E( 2 ) — 26^); these coefficients must be 
used in Eq. (|4B]) to determine the maximal occupation number of this two-boson state. The 
noncollective eigenstates have energies E( 2 ) = € N + e' N , and keep their unperturbed (har- 
monic oscillator) structure, apart from the fact that for a level Er 2 ) = 2e' N the pair wave 
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function must be orthogonal to the zero-coupled pair X^m(frjvzm) 2 - F° r 9 < O5 the two-body 
ground-state is always collective, whereas the first excited state can be either collective or 
noncollective, depending on the interaction strength. 

First we study a simple case of four bosons in four HO levels (N = 0, ... ,3), as the 
dimensionality is still sufficiently small to allow comparison with exact diagonalization. In 
Fig. [I] and Fig. ||] results are shown for attractive pairing. The weak-coupling regime is 
displayed in Fig. [I], where the exact energy is compared with the Hall-Post lower bound 
Co = 6E (2) and the new lower bound L\ = S^)Eo(2) + (6 — S^^)Ei^)- We also plot, as an 
example of a simple upper bound, the energy of the Hartree solution, which is known 
exactly for this system. If the sequence of structure coefficients wn is decreasing with N, 
then the Hartree energy E? A ^ for general A can be shown to equal E^/A = €q + (A — l)g, 
if g < g c , and E^/A = e + (A — l)g c (2 — y), if g < g c , where the critical strength is 

_ wo(ei-Eo) 
Me 2(i4.-l)(lDo+tOl)" 

Both lower bounds coincide with the exact result as g — > — 0. For more negative values of 
g, the Hall-Post bound quickly diverges from the exact result, whereas the improved lower 
bound follows the exact result quite closely. In fact, it is easy to see that the improved 
lower bound, in contrast to the Hall-Post one, becomes exact also in the strong-coupling 
limit, for all attractive pairing forces. This is because, as g —>■ —00, the interaction term 
increasingly dominates over the external potential in the ground-state energy, and the lower 
bound C\ becomes exact for a separable hamiltonian P^P. Fig. ^|, where we show the 
relative error with respect to the exact energy, demonstrates this explicitely. The kink in C\ 
at g = — 5/16 ~ 0.3 occurs because at this value of the coupling strength the first excited 
state of the A = 2 system changes from a solution of Eq.(|68|) to the unpaired solution e' + e' 1 . 

We checked that in all cases £1 is the optimal bound, i.e. S^) + $1(2) > 6> by calculating 
S™) through Eq.(^) for the second lowest solution of of Eq. (j68|) , or, if the first excited state 

is the triplet Cq 00 c\ 1im |), by realizing that the maximal joint occupation number of this triplet 
is equal to four in four-boson space. 

The results for repulsive pairing (g > 0), shown in Fig. |3], are less impressive. In this case 
we lose the feature that as g — > +00 the two-body force dominates the ground-state energy; 
the system will simply tend to make pairs orthogonal to P\ and the one-body part of the 
hamiltonian can never be neglected. Of course, the new bound C\ (which is the optimal C/x) 
is still better than the conventional Cq. 

Systems with a larger number of particles and/or shells can be similarly treated. Results 
for 1000 bosons in 50 harmonic oscillator main shells are shown in Fig. ^, and the appreciable 
improvement of the new lower bound over the conventional one is again clear. For g < g', 
where g' ~ —2.6 x 10~ 4 , we have £2 > C\ and the optimal lower bound C/x is given by 
£2 instead of C\, the difference between £2 and £± is marginal, however, certainly when 
compared with £ . 

B. Bosons interacting with power-law potentials 

As an example of a self-bound system we consider A spinless bosons in three dimensions, 
interacting with power-law potentials, 
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^=-^EV'+9 sgn(P) £ I - r, f . (69) 

i ii<i2 

Note that physically relevant potentials must have (3 > —2 to ensure an eigenvalue spectrum 
bounded from below. 

In accordance with Eq.(^), we try to derive lower bounds for the three-particle ground- 

? 0(3)[|] 



state energy _E (3)[4]; in terms of the solutions of the following relative pair hamiltonian, 



H' 2 = (~)(--V 2 )+gsgn(f3y. (70) 



.A 3/ V m 

Scaling laws can be used to write the eigenenergies E($) of H' 2 as 





3+2 



E ^ =9 {-jL-g) ^ (71) 

in terms of the eigenenergies 77(2) of —V 2 + sgn( / 9)r /3 . 

For the latter hamiltonian we determined numerically the L w = + ground-state energy 
770(2) and wave function ^0(2) — g( r )^oo, as well as the energy 771(2) of the first excited 
(symmetric) state, which has L w = + for (3 < 2 and L w = 2 + for (3 > 2. 

We also determined the maximal occupation S^f of the ground-state pair \&o(2) in + 
three-boson space, which according to Eq.(|62|) equals the largest eigenvalue A of 

(A - l)G(r) = 2 J dr'r' 2 W{r, r')G(r'). (72) 

The operator W reads as 

W(r, r') = \ £ dxg(\ r - + r' \)g(\ ± + r |), (73) 

where x denotes the cosine of the angle between r and r. We solved the radial eigenvalue 
equation (|73|) on a grid. 

The lower bounds ( |3"TD for the three-boson energy £7 (3)[4] can now be used, according 
to Eq. (|63"D , to derive lower bounds C < Ci < -Eq(A) for the general A-boson system, 



3770(2) = 3Ct7o(2) 



\Amg) 

d = C (^ 2 > {2) + [3 - S^] Vm ) . (74) 

Note that, because of the scaling properties of power-law potentials, the A-dependence of 
Co and C\ can be absorbed in the coefficient C appearing in Eq. (|74|) . 

In order to compare the new lower bound C\ with the conventional Co we have plotted 
in Fig. |5| the relative improvement R = (Ci — Co)/ \ C | , for a range of powers — 2 < (3 < 10. 
Being a ratio, R is independent of A. It can be rewritten as a product 

R = ^± = (1 _ I SS ) f "'f -*<* ) = R,R 2 , (75) 

\ Co \ \ 3 / \ I 770(2) I / 
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where the contributing factors R\ and R2 are related to the maximal occupation of the 
ground-state pair, and to the energy difference between ground and first excited state, 
respectively. These factors are also plotted in Fig. |5|. 

As can be seen from Fig. [|, the relative improvement R becomes zero for two values, 
(3 = and (3 = 2. For (3 = 2 it is R\ that vanishes, since the operator in Eq.(|7^) has 
an eigenvalue equal to one if g(r) ~ exp (— r 2 ). This reflects the fact that the conventional 
bound Co becomes exact for harmonic oscillator systems |]J. For (3 = it is R2 that vanishes, 
since the pair energy spectrum becomes degenerate as (3 — > ±0, that is, 77^2) — > ±1 in this 
limit (see, e.g., [fffi]). 

Except for extreme values of (3 (which means (3 close to —2 or positive and large) the 
improvement of C\ over Co seems modest, e.g. for the case of gravitating bosons {(3 = — 1) 
we find R ~ 1.4%. However, for most power-law potentials the conventional bound Co is 
already quite a good approximation to the exact energy of the three-body system, so any 
improvement cannot be large on this scale. In Table 1 we compare, for the three-body system 
with m = g = 1, the exact ground-state energy (taken from jlj) with the lower bounds Co 
and C\, for a few values of (3. It is seen that in the three-body system the improved bound 
removes a sizeable fraction (between 25% and 75%) of the remaining discrepancy between 
the exact energy and the lower-bound of the Hall-Post type. 

It is interesting to note p3f that in the limit (3 — > the power-law potential is related 



to the case of the logarithmic potential, lnr = lim / 3^ ( r — As a consequence, the 

eigenvalues r] v ((3) of the hamiltonian —V 2 + sgn(/3)r^ are, for small (3, connected with the 
eigenvalues rfc of the hamiltonian —V 2 + lnr via the relation 



tjM - sgn(/3)+ I (3 



Vv 



-mi 
2 







(76) 



which holds up to terms decreasing faster than linear in | /3 |. Eq.(|76|) shows the origin of 
the degeneracy in the power-law eigenvalue spectrum for (3 — > 0. This degeneracy is absent 
if we consider directly the logarithmic potential, i.e. an A-boson system with hamiltonian 



zm i h<i 2 



(77) 



A straightforward analysis then leads to 



C, 



A(A 



6 



-9 



A(A-l) 
C 1 = ^^g 



0(2) 



2 \Amg J 

max„L 1 To omaxl„L 



S mVo(2) + t 3 - S oi2)]Vi(2) + 2 ln {j^gj 



(78) 



where ^ (2) = 1.0444332, »jf (2) = 1.847442, and Sfiffi = 2.986419. To compare this with the 
values in Fig. || one can consider, e.g., mAg/2 ~ 1. The relative improvement then becomes 
R = w 0.35%. 

\Co\ 
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C. Electrons confined in a harmonic oscillator well. 



As an example of a fermion problem we improve the bounds derived recently by Juillet 
et al. for a quantum dot system of electrons confined in a harmonic oscillator well ||13| . The 
hamiltonian in atomic units (m = e = 1) reads 

i=i V z z / i<j I Li l j I 

The harmonic center-of-mass motion can be split off and treated exactly, and we concentrate 
on the relative motion. 

We consider a system of three electrons with total spin S = |, e.g. three spin-up electrons. 
The spatial wave function is antisymmetric, so the bound derived in Section [V] can be applied 
for the energy of the ground state, which has L n = 1 + . 

According to Eq. ([□]), with N = 2 and A = 3, we must first construct the ground state 
of the relative two-body hamiltonian, 

* = I (" v2 + T r2 ) + (80) 

In the previous examples the pair ground state was nondegenerate. In the present case the 
lowest antisymmetric eigenstate g(r) forms a L n = 1~ triplet, 

^(r) = g(r)Y lfl (Q). (81) 



Since the pair ground state is now degenerate, we must generalize Eq.(^9j), and maximize the 
joint occupancy A = X)„ S(2)n of the members of the triplet. This leads in a straightforward 
fashion to an eigenvalue equation, 



(A-l)G^(r) = -2^ Jdr' 



G M v(r'), (82) 



which replaces Eq.(p2]). The overlap function G^r) between the 1 pair state and one of 
the members ^(3)u of the A = 3 ground-state 1 + triplet, has the following tensor structure, 

G lu ,(r) = G(r)(l(ilu- f i\ li/)Y w (fi). (83) 

Substitution into Eq.fl8"2[) leads, after some angular momentum algebra, to a radial eigenvalue 
equation 

(A - l)G(r) = 2 J dr'r' 2 W{r, r')G(r'). (84) 

The operator W reads 

W( rr ) = _ rr / da; (P (x) - P 2 (x)) , , , r , , , (85) 

2 J-i | § + r' 1 1 | + r | 

where P; are the Legendre polynomials. 
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We checked that for a HO p wave function, g(r) ~ rexp (— r ), the maximal eigenvalue 
of Eq.(84) yields A = 3, which equals the number of pairs in the A = 3 system. This means 
that the Hall-Post bound Cq already coincides with the exact result for HO systems, as 
could be expected from the discussion in [13| . 

In the presence of the electron-electron repulsion the pair wave function is distorted from 
the HO shape, and we find a maximal eigenvalue A = 2.901644 for uj = 0.01. The lowest 
energies of the hamiltonian (BDj) are eo = 0.05560837 for the (p-wave) ground state and 
ei = 0.06171271 for the (/-wave) first excited antisymmetric state. For the Hall-Post bound 
we thus find £ = 3e = 0.1668, in agreement with (T3|. This is already quite a good bound 



compared with the exact three-body ground-state energy = 0.1680, as quoted in p3 
The new bound improves this to C± = Aeo + (3 — A)ei = 0.1674 . 

For uj = 10, which is closer to a pure harmonic oscillator system, we find A = 2.999780, 
e = 18.32273, and e x = 31.14720, yielding bounds £ = 54.968 and C x = 54.971, to be 
compared with the exact result = 54.973, quoted in [ |T3]j . 

In conclusion, by taking the structure of the pair wave function into account we are able 
to halve the remaining deviation between the Hall-Post type lower bound and the exact 
result. 



VII. SUMMARY 

Motivated by the renewed interest in lower bounds for the ground-state energy of many- 
body systems, we have developed a method to improve the existing lower bounds of the Hall- 
Post type. The method is based on an exact sumrule for the energy in terms of two-body 
occupation numbers (or, equivalently, spectroscopic factors related to A — 2-particle removal) 
in the A-particle ground state. The pair occupation numbers that enter the sumrule refer to 
the two-body eigenstates of the two-body cluster hamiltonian in the conventional Hall-Post 
decomposition of the many-body hamiltonian. We find that it is possible to derive upper 
bounds for these pair occupation numbers, without detailed knowledge of the structure of 
the A-particle wave function. These upper bounds, or maximal pair occupancies, do depend 
on the structure of the pair state, and can be used to obtain strict lower bounds to the 
A-particle energy which are better than the conventional one. 

We have studied both the bosonic and fermionic sector, and developed a framework for 
both fixed-center systems and self-bound systems, where the wave functions are transla- 
tionally invariant. We have applied the formal results to various numerical examples, and 
demonstrated that significant improvements are obtained over the conventional lower bound. 

Several problems are still remaining. In the case of self-bound systems a method to eval- 
uate the maximal pair occupation in A-particle space is not available for A > 4. Also states 
of mixed spatial symmetry are not yet treated. If the two-body energy spectrum contains 
a continuum part, the associated pair strength does not contribute to the lower bound in 
the present work; a better treatment of the continuum part would be very interesting, as 
it would lead to the derivation of improved bounds for the critical coupling strength 0-f|] 
needed to achieve binding in many-body systems. 

In general, further improvements could be made by more refined approximations to the 
distribution S u ^) of the pair occupations over the various pair states. In the present work 
this occupation is maximized for each pair state separately; in reality, of course, they are 
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interrelated, as they reflect occupations within the same A-particle state. Such a refinement 
is in particular needed for the fermion case, since the noninteracting limit is at present not 
reproduced. 

The present work can also be rephrased in terms of abstract many-body theory. The 
A-particle ground-state energy in Eqs.( p3| , p6|) is expressed as 

E 0[A) = Trace {^ 2 [a]p {2) }, (86) 

where a = 1/(A — 1) or a = 2/ A for fixed-center and self-bound systems, respectively. 
Minimizing the r.h.s of Eq. (|8~6D over all A-representable two-body densities would yield the 
exact A-particle energy. The full set of exact conditions for A-representability are of course 
unknown. Minimizing the r.h.s. of Eq.(|86"D over all two-body densities which comply with a 
limited set of A-represent ability conditions will then yield a lower bound for the A-particle 
energy. The conventional bound Co can be seen as the lowest-order approximation in this 
scheme, since only the normalization condition Trace {p^} = A(A — l)/2 is required for 
the two-body density matrix in Eq.([8^). The improved bounds in this work can be viewed 
as imposing additional conditions on the natural pair occupation numbers of the two-body 
density matrix in this scheme. 

This work was supported by the Fund for Scientific Research-Flanders (FWO- 
Vlaanderen) and the Research Council of Ghent University. 

APPENDIX A: NATURAL ORBITAL REPRESENTATION FOR PAIR STATES 

1. Two-boson states 

In a general SP basis, a two-boson state can be expanded as 

l*(2)> = 4z<^4i (ai) 

V2 aj3 

where C a p = Cp a is a (complex) symmetric matrix, and TracejCC"^} = J2 a p I | 2= 1 
fixes the normalization. 

The hermitian one-body density matrix reads 

Pap = <* ( 2 ) I c\c a | * (2) ) = 2 ]T C aX C; x , (A2) 

A 

or, in matrix notation, p = 2CC^ . 

Under a unitary transformation c^, = J2a U a a'C a to a new SP basis, the matrices p and 
C transform as 

p> = u^pU, C = U ] CU*. (A3) 

We can always make a unitary transformation to the natural SP basis that diagonalizes 
the one-body density matrix p. In this basis p is real, and as a consequence the commutators 
[C, C*] = [p, C] — vanish. It follows that in the natural basis C is block diagonal, each 
n x n block C^' corresponding to a ra-fold degenerate p aa . 
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For such a block the matrix = C^J 2/ p aa is a unitary and symmetric matrix, which 
can be diagonalized by a real orthogonal transformation. Indeed, if X is an eigenvector 
of L> (a) with eigenvalue A, then D^X = XX implies D^X* = XX*, because £> (a) * = 
and A* = A -1 . So either X = X* is real, or A is degenerate with eigenvectors 
X, X*, which can be replaced by an orthogonal pair from the real linear combinations 
X + = (X + X*), X_ = (X- X*)/i. 

Since the transformation to the basis of real eigenvectors of C is real orthogonal, it also 
corresponds, according to Eq. (|A3|) , to an allowed unitary transformation on the SP basis. 
In this basis, = 6 mn \J p aa / 2 exp (i9 n ) . The phase can be absorbed in the SP states. As 
a result, the desired canonical form of a two-boson state reads, 

l*( 2 )> = 4fE^(4) 2 l>, ( A4 ) 



where x a = Jp aa /2 is real and positive. 



2. Two-fermion states 



In a general SP basis, a two-fermion state can be expanded as 

l^(2)) = ^E^ C Wl)' ( A5 ) 

where C a p = Cp a is a (complex) antisymmetric matrix, and TracejCC^} = 1. 

An identical analysis as in the bosonic case leads in the natural basis to matrices 
which are now unitary and antisymmetric, and can be brought to a canonical form by a 
real orthogonal transformation. If D^X = XX, then D^X* = -XX* (because D^* = 
—D^ ). So the eigenvalues/vectors come in pairs (X,X), (—X,X*), where X^X* = 0. For 
such a pair we may transform from the eigenvector basis X, X* to the real basis X + = 
(X + X*)/V2,X_ = (X - X*)/iy/2. The 2 x 2 diagonal block in the X,X* representation 
is transformed as 

(A6) 

The transformation to the real basisvectors X + ,X_ is again real orthogonal, and cor- 
responds to an allowed unitary transformation on the SP basis. In this basis, = 




fimny Paa /2 exp (i6 m ), where n, n are associated pair states and 9 n = n + tt. Apart from an 
overall sign, the phase can be absorbed in the SP states. As a result, the desired canonical 
form of a two-fermion state reads, 

I *( 2) > = 4= E - 44) I) = V2 E *«44 1>. ( A7 ) 

V ^ «>0 a>0 



where x a = yp aa /2 is real and positive, and the summation a > is made over distinct 
pairs. 
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APPENDIX B: MAXIMAL PAIR OCCUPANCY FOR A GENERAL NUMBER 

OF BOSONS 



We follow here closely the reasoning by Pan et al. [20] for the fermion pairing problem. 
A two-boson state ^(2) is expressed in its natural basis as 



*(2)> = ^EM4) 2 D = W 2) I>, 



(Bl) 



where the x a are real and positive, and ^Z a x a = 1. The construction of the eigenvalues of 
the V 9 (2)V 5 (2) operator then proceeds as follows. 

The uncorrelated A-boson states in the natural SP basis read 



^)> = n(4ri> 



(B2) 



i=i 



where J2iPi — A, and can be classified according to the broken pairs they contain, i.e. with 
Pi = 2rrii + rj and = or = 1, we can construct the corresponding vacuum | 0) of the 
(f(2) operator, 



o> = nw) p< D- 



(B3) 



i=l 



Obviously, the ^(2)^(2) operator does not connect A-boson states belonging to a different 
vacuum | 0). For our purpose we may take | 0) =|) (the zero-boson state) if A = 2a is even, 
and I 0) = ct |) =| /3) if A = 2a + 1 is odd, where (3 is one of the single-boson states. 



Next we introduce generalized operators 



V2 



N(y) =E T 



X' 



c 1 c -I — 



which obey the commutator algebra 



,y\N{yi) - y 2 N(y 2 



y\ - V2 



yi - y-i 

One can now show that for a suitable choice of variables y%, . . . , y a -i, the vector 

\* (A) ) = <p\0)<pi(y 1 )...<p\y a - 1 )\0), 

is an eigenvector of the ^ (0)^(0) = ^[2)^(2) operator. 
Using the commutation relations (p5|) , one finds 



a-l 



1 



^(0)^(0) I V)> = 2A(0) +4^- I V)> 

+2E My* 

k=l 



V 



2 

Vk 



a— 1 c 



a-l 



(B4) 



(B5) 



(B6) 



(B7) 

^(0)^(0) J] ^(2/0 I 0), 



i = 1 

(i ± k) 
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where A(y) is the eigenvalue for N(y) acting on the vacuum | 0), 

A(y) = E (I + s6 ^) ■ (B8) 

For even A, | 0) =|) and s = 0; for odd A, | 0) =| (3) and s = 1. 

The second term in Eq. ([B7|) vanishes if the variables (yi, . . . , y a -i) are solutions of the 
set (|46|) of nonlinear equations, whereas the coefficient in front of the first term yields the 
corresponding eigenvalue (S3) of the <^t(0)<^(0) operator. 
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TABLES 



TABLE I. Energies and bounds in a three-body system interacting with power-law potentials 
sgn(/3)r*. The hamiltonian is given by Eq,(|69D with A = 3 and m = g = 1. We compare the lower 
bounds £q and C\ with the exact energy E / 3 \, taken from [jjj and properly rescaled. The last 
column contains the ratio (C\ — £q)/(-E'o(3) — A)). 



p 


^0(3) 


Co 


A 


ratio 


-1.0 


-1.0670 


-1.1250 


-1.1095 


27 


-0.5 


-1.4911 


-1.5043 


-1.4987 


42 


0.1 


3.6383 


3.6363 


3.6374 


57 


0.5 


5.0780 


5.0718 


5.0757 


64 


1.0 


6.1323 


6.1276 


6.1309 


71 


2.0 


7.3485 


7.3485 


7.3485 




3.0 


8.1228 


8.1163 


8.1212 


75 
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FIGURES 
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FIG. 1. Ground-state energy of four bosons with hamiltonian (^), in four harmonic oscillator 
levels. Solid line: exact result £\j(4) • Dashed line: Hartree upper bound E^. Dotted line: 
Hall- Post lower bound Cq. Dot-dashed line: improved lower bound C\. 
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FIG. 2. Relative errors for the ground-state energy of four bosons with hamiltonian 
four harmonic oscillator levels. Dashed line: Hartree result (E^ — £J (4))/-E'0(4) [where i?o(4) is the 
exact energy]. Dotted line: Hall- Post result (Co — £?o(4) )/-^o(4) • Dot-dashed line: improved result 

(A - -E' (4))/£ , (4). 
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FIG. 4. Ground-state energy -Eo(A) f° r A = 1000 bosons with hamiltonian (U) in 50 harmonic 
oscillator levels. Dashed line: Hartree upper bound. Dotted line: Hall-Post lower bound Cq. 
Dot-dashed line: improved lower bound C^, where jx = 1 or fl = 2 (see text). 
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FIG. 5. Lower bounds for the energy of a system of A bosons interacting with a power-law 
potential sgn(/3)r /3 . Solid line: relative improvement R = (C\ — Co)/ \ Co I- The ratio R = R1R2 
is a product of two factors (see Eq.([75|)), which are also plotted. Dashed line: R\ = 1 — 5^2*/3. 
Dotted line: R 2 = {t] 1{2 ) - %(2))/ I %(2) I- 
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